function imP = polar (I)

[F,C] = size(I);
Om = (F+1)/2;           
On = (C+1)/2;
sx = (F-1)/2;           
sy = (C-1)/2;
img  = zeros(F,C);
delR = 1/(F-1);
delT = 2*pi/C;
rmin=50;
rmax=150;

for ri = 1:F
    for ti = 1:C
 
       r = (ri - 1)*delR;
       if(r>0.2 && r<0.4)
         t = (ti - 1)*delT;
         x = r*cos(t);
         y = r*sin(t);
         xR = x*sx + Om;
         yR = y*sy + On;
         %img (ri, ti) = I(round(xR), round(yR));
         img(round(xR), round(yR))=1;
       end
    end
end
imP=img;
end


